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Abstract 

We present spectra of accretion discs around white-dwarfs calculated with an improved and updated version of the 
Shaviv & Wehrse (1991) model. The new version includes line opacities and convective energy transport and can be 
used to calculate spectra of hot discs in bright systems (nova-like variables or dwarf novae in outburst) as well as 
spectra of cold accretion discs in quiescent dwarf novae. 
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1. Introduction 

In weakly magnetized cataclysmic variable stars 
(CVs) matter lost by a Roche-lobe filling low-mass 
companion forms an accretion disc around the 
white-dwarf primary. In bright systems such as 
nova-like binaries and dwarf-novae in outburst the 
disc is the dominant source of luminosity but even in 
quiescence the disc emission might provide crucial 
information about the physics of accretion, in par- 
ticular about the mechanisms transporting angular 
momentum and releasing gravitational energy. 
Until now one could c alculate only spec t ra for 



hot accretion discs (see IWade fc Hubenv , 1998 



and references therein). In such an approach the 
disc is divided into rings whose vertical str ucture is 
calculated u sing the program TLUSDISK Irlubenv 
(|199(A Il99ll) . A different way of solving radiative 
tr ansfer equations in accret ion discs was presented 
bv lShaviv fc Wehrse I (|l99ll hereafter SW). Also, in 
their approach one could apply it to hot discs only. 
The preliminary attempt to apply the SW code to 



cold quiescent discs in Idan et al. ( 1999h was not 
conclusive. The main obstacle in solving the verti- 
cal structure of cold accretion discs is the presence 
of convection-dominated zones . (The description of 
the quiescent disc of SS Cyg bv lKromer et al. , 2007 , 
neglects convection). 

The absence of spectral models for cold quiescent 
discs is the reason why fitting models to observations 
can be a frustrating exercise(see lUrban fc Sion . 
20061 . and references therein). Using the model s of 



hot stationary discs of Wade fc Hubenv ( 1998h to 
describe the spectra of cold (and non-stationary) 
quiescent discs cannot be expected to give relevant 
results. In fact since the observational data were the 
UV spectra obtained by the IUE, it was not very 
surprising that the best fits were obtained with no 
disc emission at all: quiescent discs are not supposed 
to be UV emitters so that the only source of UV ra- 
diation in such systems should be the white dwarf. 

The description of the quiescent state of the 
dwarf-nova cycle is the Achilles heel of the widely 
accepted disc instability model (DIM) according to 
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dwarf-nova outbursts are due to a thermal-viscous 



instability (see lLasotal 120011 for a review). For ex- 
ample while the DIM predicts the quiescent disc to 
be optically thick (especially in its outer regions), 
the observations of the quiesc ent dwarf-nova IP Peg 



seem to suggest the opposite (jLittlefair et all . 12001 



Ribeiro et all 120071 ). Correct modeling of the emis- 



sion from quiescent discs is therefore an obvious 
prerequisite to solving this thorny problem . 

In general a full description of dwarf-nova out- 
burst cycle should include the radiation transfer 
equations. In practice this much too costly and 
therefore unfeasible. However, instead of includ- 
ing one can combine the DIM with the solutions 
of the radiative tran sfer equations f o r the same 
disc parameters. The lHameurv et al.l (| 19981 . here- 
after HMDLH) version of the DIM in which the 
time-dependent radial evolution equations use as 
input a pre-calculated grid of hydrostatic vertical 
structures (a 1+1D scheme) is well suited to such 
cntcrprize. The vertical structures are calculated 
using the standard equations of stellar structure or 
the grey-atmosphere approximation in the optically 
thin case. The angular-momentum viscosity trans- 
port mechanism is d escrib ed by the a-ansatz of 



Shakura &: Sunvaev I (|1973[ K For a given a, M/R 



and effective temperature T e g there exist a unique 
solution describing the disc vertical structure. Such 
solutions which are very handy for solving the time- 
dependent equations of the disc evolution cannot 
be used to produce emission spectra. These can be, 
however, calculated by using the same input pa- 
rameters (a, M/R 3 and T e g) in a radiative-transfer 
code on the condition that the vertical solutions 
calculated by the two methods are the same. In this 
way one can reproduce spectra of the whole cycle 
of a dwarf-nova outburst. 

The outline of the present article is as follows. 
In section [5] we discuss the model used with spe- 
cial stress put on marking the differences with orig- 
inal SW code on which it has been based. In section 
[3] we present solutions obtained for various types 
of physical set-ups. In subsection 13.21 we compare 
our solutions for vertical structure with those ob- 
tained with HMDLH code. Then in 13.51 we corn- 
par e^ur_djsc_^p_ectra with t he UV spectra obtained 
bv IWade &: Hubenv I (|l998h . In the following sub- 
section subsec:non-stat we present spectra of cold 
non-equilibrium discs which in the DIM represent 
the quiescent state of the dwarf- nova outburst cycle. 
Section |4] contains the discussion, and conclusion. 



2. Assumptions, equations and methods of 
solution 

The structure of the cod e is based on the model 
of Shaviv fc Wehrse ( 199ll ) . The principal improve- 
ments consist in using (modern) line opacities and 
in a different treatment of the convective energy 
transport. It seems that our code is the only work- 
ing accretion-disc radiative-transfer code which 
includes convection self-consistently and is able to 
handle dominant convective fluxes. 

2.1. Assumptions 

The accretion disc is described as composed of 
concentric rings orbiting the central gravitating 
body of mass M with the Keplerian angular fre- 
quency 



n K = 



GM 



(1) 



We assume that the disc is geometrically thin, i.e. 



z < R 



(2) 



where zq is the disc height. This assumption allows 
to decompose the disc equations into their vertical 
and radial components. In this paper we will solve 
only the vertical accretion-disc equations and ob- 
tain the corresponding emission spectra. However, 
since the radiation of a steady-state accretion disc 
can be considered as a sum of radiation from indi- 
vidual rings our results can directly applied to such 
configurations. The same is true of quasi-stationary 
phases of dwarf-nova outbursts. Also radiation from 
non-steady quiescent discs of dwarf-nova stars can 
be treated as being the sum of the emission of indi- 
vidual rings. This is true as long as radial gradients 
are much smaller than the vertical ones. During the 
outburst of dwarf nova when the temperature and 
densities fronts are propagating through the disc this 
assumption might no longer be valid. 

2.2. Equations 

The vertical structure of a ring at a distance R 
from the center is found by solving: 
- The hydrostatic equilibrium equation 



dP 

dz 



-P9z 



-p£l K z, 



(3) 
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The energy balance equation 

= q = (¥k\ + (**± 

dz VIb \ dz / rad \ dz 



(4) 



where F z is the energy flux in the vertical 
(z) direction. We u se th e a prescription of 
Shakura fc Sunvaev I (|19731 ) to describe the vis- 
cous energy generation 



3 o P 
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Q vis (z,R) 
The radiative energy flux is given by 



(JO, R, A) - B(T(z, R), A)) K (\)d\ 



(5) 



(6) 



= -Qvis(Z)-E) 

where J is the mean intensity, B the Planck func- 
tion and k(A) is the monochromatic mass absorp- 
tion coefficient. 

The radiative transfer equation is 



n • VI = k(S - I), 



(7) 



where I is the intensity, S the source function and 
n is a unit vector in the ray direction. 
Following HMDLH the con yective energy flu x is 
calculated by the method of lPaczvhskil ( 19691 ). 
Whenever the radiative gradient 



Vrad = 



d\uT 
d\nP 



(8) 



rad 



is superadiabatic, the temperature gradient of the 
structure V is convective (V = V con v)- The con- 
vective gradient is calculated in the mixing length 
approximation, with the mixing length taken as 
H m i — a m \Hp, where Hp is the pressure scale 
height: 



H 



P 



(9) 



pg. + iPp) 1 ^' 
The convective gradient is found from: 

v conv = v ad + (V rad - V ad )F(y + A) (10) 

where V a d is the adiabatic gradient, and Y is the 
solution of the cubic equation: 



9 
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-Y A + VY Z + V Z Y - V = 



(11) 



' ml 



where r m i = npH m \ is the optical depth of the 
convective eddies. It is assumed that the eddies 



are optically thick so that the convection is effi- 
cient, yet the contribution to the adiabatic gra- 
dient V a d + (V ra d — Va d )y(y + A) cannot be 
neglected. In those cases in which the convection 
reached the photosphere the effect on the weak 
lines must be noticeable. The coefficient V is given 
by: 



v- 2 = 



'ml 



3r m i 



512a 2 T 6 H P 
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dWT 



(Vrad - V ad ) (12) 



In our model we take a m i = 1.5 
The outer boundary condition is 

Qvis = ""^eff 

In the case of a stationary disc 

z 



(13) 



(14) 



M = 2ttR J pv r dz = constant, 
o 

where v r is the radial flow velocity one has 



&7TO~ \ R 3 



R 



but our solutions apply to any distribution of ac- 
cretion rate M(r) in a geometrically thin accre- 
tion disc. 



2.3. Input (opacities and EOS) 

We have improved several features of the original 
SW code but the use of modern opacities and espe- 
cially of the line opacities is the major improvement 
that we present here in some detail. 

There are several ways of calculating opacities. 
They can be obtained using e ither codes such as 



( Kurucz 



PHOENI X (lHauschildt et all . Il999h . "Atlas 12" 
1992TT or from the database of the OP 



project (jSeaton et al 



199J). The updated version 
of the SW code can use the opacities from "At- 
las 12" but since calculating the opacities is time 
consuming we decided to use tabulated data from 
the OP database. To tabulate the data we used 
the data from the Opacity Project - OPCDJ5.3 
that were downloaded from the website at the Cen- 
tre des Donnes de Strasbourg (CDS). For a given 
abundance mixture the subroutines mx.f, mixv.f, 



3 



opfit.f and mixz.f allowed us (with some modifica- 
tions) getting the tabulated opacities for the wave- 
lengths chosen as our basic grid. For the emerging 
spectra from every ring we use 10,000 wavelengths 
starting from 100A up to 10 6 A. The grid interval 
varies between AA = 0.25A for the most important 
zone such as the range between 850 - 2300A and 
1 — 2A used for wavelengths shorter than 800A or 
above 2300A. For wavelengths longer than 5000A, 
AA varies logarithmicaly starting at 2A and getting 
up to 200 A for the very long wavelengths. All the 
data were tabulated using a grid for the tempera- 
ture T and electron density Ne. The indices ite and 
jne are defined by 



ite = 40a; log(T) , jne = 4x log(ATe) 
A(ite) = A (jne) = 2 



(16) 



The mesh in the tables is calculated for 140 < 
ite < 320 therefore the lowest temperature that 
can be used in the code is 3160K. Let us stress 
that the opacities in the present work do NOT 
take into account the global line broadening due 
effects such as velocitie s, (expansion opacity, see 
Shavivfc Wehrsel . l2005h . the inclination of the disc 
or microscopic such as microturbulence. Conse- 
quently, the lines provide only an indication of the 
presence of a given ion and cannot be used, so far, 
for abundance determinations etc. These broaden- 
ing^ will be included in the code in the near future. 
Therefore the results presented here should be 
treated as more qualitative than quantitative in the 
sense that for the moment we will not discuss the 
ratios between the emission or absorption lines but 
rather their presence or absence in the spectrum. 

The reason for allowing the code to work either 
with OP tabulated data or with "Atlas 12" was to 
test the tabulated opacities. In Figure[T]we show the 
comparison of opacities as function of wavelength 
calculated for T = 10 5 K and p = 2 x 10~ 8 g cm -3 , 



and T = 3 x 10 4 K and p = 2 



10 7 g cm 3 using 
either "Atlas 12" or the OP. The "Atlas 12" con- 
tinuum opacities used were for a solar mixture with 
the addition of only H and He lines. The opacities 
from the OP Project are for also for mixture of solar 
abundance (which includes the lines and the contin- 
uum) . 

The Rosseland mean opacities are taken from 
OPAL tables for solar composition. While solving 
the radiative transfer equation we checked for the 
consistency between the Rosseland mean opacity 
obtaining from the solution of the radiative transfer 
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Fig. 1. Comparision between the opacities as function of 
wavelength obtained from Atlas 12 and those obtain from the 
OP for T= and p = (upper figure) and T= and p = (lower 
figure). The Atlas 12 opacities are the continuum opacities 
with only H and He lines, while the OP opacities are for 
solar mixture 

equation and the value from the OPAL table. The 
eq uation of state is interpolated from the tables 



of iFontaine fc Graboske fc Van Horn I (|1977T ). The 



equation of state used in the OP project is differ- 
ent than the one use by Kurucz or in OPAL. We 
are aware of the this problem and its effect on the 
opacities but we will not discussed it here. 

2.4. The method of solution 

The basic idea behind the SW code was to couple 
the hydrostatic with the radiative-transfer equa- 
tions. Therefore the code iterates the two equations 
simultaneously. The iteration continues until a 
structure is found to which both the solution for 
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the hydrostatic and the radiative-transfer equations 
converge. 

The temperature profile is determined from Eq. 
([6|). The temperature T(z,R) in this equation is 
solved by expressing J in terms of B using the ra- 
diative transfer equa t ion (for detail d es cripti on see 



Wehrse I (jl98lM , As 



Kalkofen fc Wehrse I (|1984|) 
mentioned above we assumed that each radial ring 
is independent of its neighboring rings hence there is 
only a vertical radiative flux. The radiative transfer 
equation is solved in the two stream approximation. 
J + represent the specific intensity in the outward di- 
rection and I~ is the specific intensity in the inward 
direction. The radiative equation is then written as: 



± 



dI ± {R, z, X) 

dz 



-(k(A) +a{X))I ± (R,z,X) (17) 
+ a(X)J(R, z, A) + k(X)B(T(R, z), X) 



where er(A) is the scattering coefficient. We use the 
same N grid points as in the hydrostatic calcula- 
tion. The boundary conditions are: no incident flux 
{IJi = 0) and at the equatorial plane the upward 
and downward specific intensities are equal {Ijpf = 
IJn)i wnere j designates the ring number. 

When solving the vertical-structure equations of 
an accretion disc the main problem is that contrary 
to stellar atmospheres the exact height of the pho- 
tosphere is unknown a-priori. Because of gravity in- 
creasing with height a change in the position of the 
photosphere is followed by a change of the gravita- 
tional acceleration which leads to a change in the 
entire vertical structure of the disc. One has to re- 
member that the upper part of the disc is optically 
thin and may include, depending on the assump- 
tions, a non negligible energy production. However, 
since this optically thin region is a poor radiator 
and absorber a small energy production there might 
have large effect on the final temperature structure. 
Once the line opacities are taken into account the 
importance of tuning the disc height is crucial since 
otherwise the outgoing radiative flux might not be 
equal to the total energy flux of the ring. Conse- 
quently the great sensitivity of the emission lines to 
the details of the temperature at the edge of the pho- 
tosphere might be used to investigate the properties 
of the energy-release ( "viscosity" ) mechanism. 

In view of the above we approach the prob- 
lem in the following way: we guess first an initial 
temperature-optical depth relation and assume 
an initial height zq. The optical depth in the disc 
uppermost point is chosen to be 10 -6 , and this as- 



Assume a T(t) law 



Guess Zq - the height of the photosphere 



Solve the hydrostatic equ -'get P and p 



Does the model produce the right energy ? 
1 No 



Make another guess of Zq 



Solve the R.T. equ. + energy equ. -*get T and z 



Is the new T[t:) - the old T(t:) The end 



Fig. 2. The iteration scheme 

sumption is checked again once the final structure is 
obtained. With the help of the T(r) relation we can 
integrate the hydrostatic equation from zo down to 
z = 0. One should emphasize that r in the current 
work is the Rosseland mean and not the optical 
depth at 5000A as in SW. When the height z is 
the correct height, the total energy generation rate 
in the ring calculated from the energy equation is 
equal to the total luminosity of that particular ring. 
If this is not the case we alter zo keeping the T(r) 
relation fixed, until the this condition is satisfied. 
Once this iteration is solved the hydrostatic model 
is consistent with the energy generation but not yet 
with the radiation field. 

We then turn to solving the radiative transfer 
equation using the disc height, the pressure and the 
density structure obtained from the hydrostatic it- 
eration. This structure is kept fixed during the it- 
eration for the radiative field. We now solve the ra- 
diative transfer equation and obtain a new T(r) re- 
lation, where new values of the temperature are as- 
signed to old values of the optical depth. If the new 
T(t) relation from the radiative iteration agrees (to 
a chosen degree of accuracy) with the old T{t) rela- 
tion from the hydrostatic iteration the hydrostatic 
structure is consistent with the radiation field. If not 
we iterate again with the new T(r) relation. 

To summarize we have 3 basic iteration loops (see 
figure [2|) in the code : 

(i) Calculating the hydrostatic structure using a 
given T(r) relation. The output results are the 
disc height, pressure and density structure. 

(ii) Solving the radiative transfer for a given disc 
height, P and p. The output result is a new 
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T(t) relation, 
(iii) Repeating the previous two iterations until 

the T(t) relation from the hydrostatic part is 

equal to the T(t) relation from the radiative 

transfer equation. 
Only after these iterations converge the emerging 
flux from the disc is calculated using the radiative 
transfer equation. 

As mentioned before the entire disc is divided into 
a series of concentric rings and their width is deter- 
mined according to their distance from the WD. All 
rings which are further away than 2i? wc i are assumed 
to have a width of li? wc i, below 2i? w d the width of 
the ring is taken to be 0.05i? w d, which allows a de- 
tailed calculation of the emission of the rings near 
the boundary layer. 

The vertical structure of the ring itself is divided 
into 100 z points (100 layers) and we use 10000 wave- 
lengths. The above procedure is repeated for each 
ring. The emerging fluxes from the rings are then in- 
tegrated in order to get the total flux from the entire 
disc according to: 



F, 



total 



F(R)RdR 



(18) 



Ru 



3. Results 

3.1. Comparison with HMDLH vertical structures 

As we mention in the introduction a 1+1D scheme 
was developed by HMDLH to describe the time- 
dependent evolution of the dwarf nova outbursts. 
The scheme uses as input a pre-calculated grid of 
hydrostatic vertical structures. For a given a, M / R 3 
and effective temperature T e g there exist a unique 
solution describing the vertical structure of a disc in 
thermal equilibrium. HMDLH code cannot be used 
to produce emission spectra. These can be, however, 
calculated by using the same input parameters (a, 
M / i? 3 and T e g ) in the present code (hereafter ILHS) 
on the condition that the vertical solution calculated 
by the two methods are the same. In this way one 
can reproduce spectra of almost the whole cycle of 
a dwarf-nova outburst. This procedure makes sense 
only if the structures calculated by the two methods 
are practically the same. 

In order to compare the solutions obtained by 
HMDLH with those calculated with ILHS we show 
in Figures [3] and [4] the vertical structures for hot 




Disc Height [cm] 



Fig. 3. Comparison between the temperature and density 
profile at R=5i? w( j and T e fj=15000K calculated by the 
HMDLH disc model and the present disc model taking into 
account the radiative transfer equation for a = 0.3 and 
a = 0.03. 

and cold discs with effective temperatures 15000K 
and 5000K respectively. The mass of the white dwarf 
is 0.6Mq, the ring radius is 5i? w d wher e the white 
dwarf radius was calculated according to lNauenberg 
(|1972r ) mass-radius relation. We used two values of 
the viscosity parameter a = 0.03 and 0.3. One sees 
a perfect agreement for the a = 0.3, while the minor 
differences for the a = 0.03 are of the order of 3%. 



3.2. S- curves 

For a given mass and distance from the center the 
thermal equilibria of accretion discs can be repre- 
sented as a T e g (E)-relation. These equilibria forms 
an S shape on the (S, T e g) plane. Figure [5] shows 
an example's of £ — T e ff S curves obtained with 
both the HMDLH code and the radiative-transfer 
calculation of the ILHS. a is constant and the ra- 
dius is 5i? w d. Each point on the S-curve represents 
a thermal-equilibrium at that radius. The upper 
branch of the S-curve is the hot stable branch, while 



6 



11000 




Disc Height [cm] 



Fig. 4. Comparison between the temperature and density 
profile at R=5iJ w( j and T c g =5000K, a = 0.03 calculated by 
the HMDLH disc model and the present disc model taking 
into account the radiative transfer equation. 



the lower branch where the temperature is below ~ 
6000K, represent the cold stable solutions. The mid- 
dle branch of the S-curve corresponds to thermally 
unstable equilibria. 

The agreement between the two codes is very good 
especially on the upper branch of the S-Curve. How- 
ever, the disc at the lower turning knee of the S-curve 
(in the unstable zone) is fully convective. Here we 
were unable to calculate the structure with the full 
radiative transfer equation code. The degree of ac- 
curacy that was chosen for the T — r to converge in 
these calculations was 1%. But since near the lower 
turning knee of the S-curve the disc is almost com- 
pletely convective and any small change in the ini- 
tial conditions such as the initial optical depth will 
affect the structure of the disc, we assume that the 
acceptable error are of the order of few percent. The 
agreement between the two codes on the lower sta- 
ble branches of the S-curve is very good. 
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Fig. 5. S — T c ff curves for R = i? w( j, Af wd = O.6M0, and 
a = 0.03 and 0.3. The solid line is the HMDLH model while 
+ represent the results from the ILHS disc calculation. 

3.3. The single ring spectrum 

In Figures O [7] and [5] we present the spectra ob- 
tained from single-ring calculations at i? wc i = 5R.0 
for M w( j = O.6M . The spectra (in arbitrary units) 
were calculated for a cold disc with a = 0.03 T e g = 
5000K and a hot disc with a = 0.3, T cS = 15000K. 
The spectrum of the cold disc is dominated by nar- 
row emission lines while the spectrum of the hot 
disc is dominated by wide absorption line. In the 
hot disc spectrum ([5]) we can identify the Balmer 
lines (the 4100A, 4340A, 48601 lines). The absorp- 
tion lines have various structures : in some cases such 
as the Lyman lines or the H y line in hot disc, there is 
an emission line in the center of the absorption line. 
We should emphasize again that one must be care- 
ful when interpreting the results at this point, since 
in this paper the opacities do not include line broad- 
ening and are not taking into account the Gaussian 
instrumental broadening function that affects and 
smears most of the fine lines. Such broadening can 
change or even eliminate some of the fine structure 
of the lines that is seen in the figures. 

3.4. Stationary disc solutions for high accretion 
rates and a — 0.3 

In Figures [9] we present an example for the total 
spectrum of a hot disc a = 0.3,M wc i = O.8M and 
accretion rate of 10 17 gs -1 . We show spectra for in- 
creasing disc area - from the inner part of the disc 
(the inner 2R W ^ area) up to the largest outer radius 
calculated - 35i? w d. One can clearly see the different 
area in the disc in which different type of emission 
is formed. For example: the Balmer jump which will 
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Fig. 6. The spectra obtained from a single-ring calcu- 
lation for cold (a = 0.3,T cff = 15000K) and hot disc 
(a = 0.03, T e ff = 5000K). The ring parameters are 
H wd = 5Rq and Af wd = 0.6M Q . 
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3500 3600 3700 3800 3900 4000 4100 4200 4300 
Wavelength [A] 



Fig. 7. The flux at the wavelength between 3500 - 4000A, 
for cold disc (a = 0.03, T cS = 5000K). The spectrum is 
dominated by emission lines. 

dominate the total flux from the disc is formed on 
the outer part of the disc, while the inner part of the 
disc is dominated only by few absorption lines. 

3.5. Comparison with Wade & Hubeny hot 
accretion disc spectra 



Wade fc Hubeny I (|l998l hereafter WH) calcu- 



lated a large grid of far- and mid- ultraviolet spectra 
(850-2000 A) of the integrated light from steady- 
state accretion disks in bright cataclysmic variables. 
All their models are for disc atmospheres which are 
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Fig. 8. The flux as function of wavelength, for hot disc 
(a = 0.3, T c g = 15000K). The spectrum is dominated by 
absorption lines and in some cases narrow emission line in 
the central of the absorption "well" 




Fig. 9. The integrated disc spectrum as function of the disc 
radii for a hot disc a = 0.3, M wt j = 0.8Mq 

optically thick and are viewed from a distance of 
lOOpc. 

Figure [10] shows a comparison between a WH ac- 
cretion disc spectrum (model cc) and a spectrum ob- 
tained with our model. The model parameters were 
are: M wd = 0.8M Q and logM = -S-hM^ 1 . In 
our disc model we used and the same number of 
rings as in their paper. The viscosity prescription 
in WH is different from the a ansatz used in our 
model. Their viscosity is based on R eynolds num- 
ber approach ( KH z fc Hubenv . 19861 ) . A value of 
Re = 5000 was assumed for all radii. We estimated 
that a = 0.3 corresponds to the parameters of WH 
model cc. 

WH spectra are smoother than our spectra since 
their models have been convolved with a Gaussian 
instrumental profile but one can expect that in our 
spectrum most of the narrow emission lines will dis- 
appear after such convolution. In general both spec- 
tra are similar. 

3.6. Quiescent dwarf-nova disc: on- stationary 
solution a = 0.03, T cS = 5000. 



Fig. 10. Comparison between Wade and Hubeny model cc 
spectrum (thick black line) and the spectrum from our disc 
model (grey line) for a = 0.3 and logAf = — 8.5M0yr — 1 . 

The total flux for an outer radius of 25i? w d is 
shown in figure QT] The spectrum between 3000 — 
10000A is nearly a flat as observed, no Balmer jump 
is present. As can be seen from figure lTTI for this par- 
ticular calculation, there is no contribution from the 
disc spectrum in the UV range. 



4. Future work 

In the future we will produce a grid of disc spec- 
tra will be used to make comparisons between ob- 
servations and the model spectra. In particular will 
focus on the quiescence disc and the effect of the ir- 
radiation by the white dwarf and the formation of a 
corona. Combining ILHS with the time dependent 
HMDLH code will provide a tool to study the dwarf 
nova outburst cycle. 

The variations in the shapes of selected lines that 
are formed in different parts of the disc will be 
studied with the purpose to develop spectral-disc- 
tomography. 



Dwarf-nova quiescent discs are non-stationary, i.e. 
M ^ constant since they are filling-up getting ready 
for the next outburst. The whole disc must be cold, 
its temperature everywhere must be lower than the 
critical instability temperature. In general, for most 
of the quiescence the temperature profile is roughly 
flat. We modeled as an example the spectrum of such 
a disc assuming T g — 5000K (which is a bit high 
for a real disc since the critical temperature is close 
to that value). 
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